Probing the largest cosmological scales with the CMB- Velocity correlation 
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Cross-correlation between the CMB and large-scale structure is a powerful probe of dark-energy 
and gravity on the largest physical scales. We introduce a novel estimator, the CMB-velocity 
correlation, that has most of his power on large scales and that, at low redshift, delivers up to factor 
of two higher signal-to-noise ratio than the recently detected CMB-dark matter density correlation 
expected from the Integrated Sachs- Wolfe effect. We propose to use a combination of peculiar 
velocities measured from supernovae type la and kinetic Sunyaev-Zeldovich cluster surveys to reveal 
this signal and forecast dark-energy constraints that can be achieved with future surveys. We stress 
that low redshift peculiar velocity measurements should be exploited with complementary deeper 
large-scale structure surveys for precision cosmology. 
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INTRODUCTION 

The Integrated Sachs- Wolfe (ISW) effect carries infor- 
mation on the evolution of the universe at low and mod- 
erate redshifts (z < 3) and thus is a powerful probe of 
dark-energy (DE) properties [1-3]. Since the ISW effect 
is most sensitive to gravity on the largest scales (hun- 
dreds of Mpc), it also offers a novel way to distinguish 
dark-energy from modified gravity theories [4-6]. Mea- 
suring the ISW signal that contributes to the total CMB 
temperature power spectrum is mainly limited by sam- 
pling variance errors on large scales, that has tradition- 
ally challenged its detection. An alternative path to ex- 
tracting the ISW effect is by cross-correlating CMB maps 
with tracers of the large-scale structure, such as galaxies. 

Recently, first detections of the ISW effect from cross- 
correlation analyses of WMAP data with different galaxy 
surveys, and x-ray maps have been obtained [7-13] . Re- 
ported detections are at the 2 — 4 a level and favor a 
DE dominated universe, independent of Supernovae type 
la (SNe hereafter). Future galaxy surveys (such as DES, 
Pan-STARRS, WFMOS) hold the promise of raising the 
significance of these detections by going deeper and wider 
to optimally sample the ISW signal. Primary anisot- 
ropics of the CMB are the dominant source of noise in 
such correlation analyses that seriously limits the ability 
with which this probe can constrain cosmology [14, 15]. 
In addition, CMB-galaxy correlations are affected by sys- 
tematics such as galaxy bias and shot noise. Moreover 
reconstructing the gravitational potential (and its time 
evolution) from the density field is a noisy process in- 
volving second order spatial derivatives. 

Alternatively, velocity flow measurements from esti- 
mates of the luminosity distance and redshift are unbi- 
ased tracers of the gravitational potential, and the latter 
can be simply reconstructed through first spatial deriva- 
tives of the previous. Also, these measurements are not 



affected by shot-noise unlike galaxy densities (specially 
in deep samples). Peculiar velocities have more power 
on the largest (linear) scales than the dark matter den- 
sity, what follows from the continuity equation (see be- 
low), and thus provides a natural counterpart to the 
density field in reconstructing the gravitational poten- 
tial. However probing large-scale velocity flows is not 
an easy endeavor. Although galaxy surveys now sam- 
ple large volumes, the intrinsic inaccuracy plaguing dis- 
tance estimations (20-25%) limits the volume available 
for such measurements [16-18]. Nonetheless better dis- 
tance indicators and other probes of the velocity fields are 
already used in cosmology but their power to reveal the 
ISW effect and the physics behind it has been neglected 
so far: SNe currently yield distance measurements with 
a 5%-10% error [19, 20] and kinetic Sunyaev-Zeldovich 
(kSZ,[21]) could in principle deliver redshift independent 
velocity errors around 100 km.s -1 . These surveys al- 
ready allow to map the large-scale velocity field at low z 
with a good agreement with other probes [19, 22, 23]. 

In this paper we propose to use the large-scale veloc- 
ity flows from low-rcdshift SNe and kSZ cluster surveys 
to measure the ISW effect from the CMB-velocity corre- 
lation. We propose to include velocity measurements at 
low redshift in combination with deeper large-scale struc- 
ture probes in tomographic analyses to better probe DE 
and gravity on the largest scales. The complementarity of 
velocities with respect to densities comes from the larger 
signal-to-noise of the CMB-velocity correlation estima- 
tor at low redshift (up to a factor of 2 gain depending 
on cosmology), as we will show in § below. Unless oth- 
erwise stated, we will focus on flat ACDM models with 
n DE = 0.75, n b = 0.05,n s = l,h = 0.7,a 8 = 0.9, and as- 
sume large-scale structure sources follow a density distri- 
bution dn/dztx z 2 exp [— (z/ (\/2z m )) 1 ' s '~\ , with a width, 
a z ~ z m /2, where z m is the survey median redshift. 
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FIG. 1: Cross power spectra for the CMB temperature- 
velocity, Tv (top), and the temperature-density, TS (bottom), 
for different mean survey depths, z m . Most of the signal in 
the Tv correlation comes from larger scales (lower multipoles) 
than the TS. We amplify the velocity correlation signal by 100 
to match the ISW amplitude from the density correlation. 



MODELING THE OBSERVABLES 

Dynamics of the fluctuations of the cosmological fields 
on the largest scales can be accurately described by lin- 
ear perturbation theory. In this regime, the Poisson and 
continuity equation in Fourier space are given by (Pee- 
bles 1980, §2) -k 2 /a 2 <$>{k,t) = 3/2H 2 nD5(k,t a ), where 
we have used that H 2 = (87rG/3)p/0, with p being the 
mean matter density, a(t) is the FRW scale factor, k is the 
comoving wavenumber and D(t) accounts for the linear 
growth of dark-matter density perturbations. The linear 
continuity equation reads, v(k,t) = —iDaS(k,to)/kk 
where k = k/fc and • denotes derivative with respect 
to comoving time t. Combining both Fourier equa- 
tions we get the linear evolution of velocity modes as 
a function of the present day gravitational potential 
{e.g. [24]), v(fc, t) = i k/H $(fc, t )g(t)k where $(fc, t) = 
D/a <b(k,to) is the time evolution of gravitational poten- 
tial modes and g(t) = 2a 2 /3Qo\/Qm/a 3 + dD/da. 

Angular cross-power spectra between the CMB tem- 
perature anisotropy and the dark-matter density S and 
peculiar velocity v fields can be derived from the cor- 
responding angular 2-point correlation function wxy = 
(A(f)F(f')) = £(2^ + l)/4rr Cf Y P e (r • f ), and yields, 



nXY 



^ J k 2 dk Pw(k)W e x (k)Wf{k). 



where the 3D 
(2n) 3 V X Y(k)5(k 



power spectra is defined 
k') = (A(k)f(k')), and 
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field spreads in comoving scales are given by ([25]) 
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kernels W*(k) that define how the signal from the X 



being dn/dr] the space density of sources and we denote 
the radial velocity v r — |v(r) ■ r\ = v. On sufficiently 
large (linear) scales, gravitational instability generates 
irrotational flows, thus we expect vanishing tangential 
component v± = 0, (only the line-of-sight modes k = kn 
are non-zero, [26]) and we define j'(x) — dj(x)/dx with 
je(x) being the spherical Bessel function of order I. In 
our computations, we shall use velocities in dimensionless 
units, by normalizing to the speed of light. 

Fig. 1 shows the CMB temperature- velocity, Tv, power 
is about 2 orders of magnitude smaller than that of the 
temperature-density, TS, estimator. However, as we shall 
see below, the noise contribution is more than 100 times 
smaller than for the TS, what makes this signal easier to 
extract in cross-correlations analyses. Note that most of 
the signal in the Tv correlation comes from larger scales 
(lower multipoles) than the TS and therefore very wide 
surveys are required to measure the expected correlation. 
On the other hand, as illustrated by Fig. 2, the peak 
amplitude of the Tv kernels (i.e, integrand in eq. (1)), 
for relatively shallow surveys, z < 0.5, peaks at I = 2, 
whereas for deeper surveys the signal is dominated by 
higher £'s. Note that the CMB-velocity is mainly sensi- 
tive to modes at k ~ 0.01 — 0.001 Mpc~ , i.e, almost one 
order of magnitude larger comoving scales than for TS. 

OBSERVATIONAL PROSPECTS 

Current probes of the velocity flow on large scales in- 
clude galaxy peculiar velocities, SNe galaxy host pecu- 
liar velocities and the kinetic Sunyaev-Zel'dovich effect 
(kSZ). The radial component of peculiar velocities of 
galaxies and SNe, v is determined by inferring their dis- 
tance, d, and then subtracting off the Hubble flow con- 
tribution to the measured redshift, z, as v = cz — Hod 
[34]. The empirical correlation -between light-curve 
shape and luminosity, and between color-curve shape and 
extinction- used to infer the SNe luminosity distance cur- 
rently yield a dispersion in apparent magnitude, m, of 
<r(m) = 0.1 — 0.15. Since m is related to the luminosity 
distance in Mpc, dh, and the absolute magnitude, M, as 
m =51og 10 d L + M, a dispersion of a(m) = 0.1 entails 
a redshift independent distance measurements errors of 
5% when ignoring the effect of the marginalization over 
the constant M. Therefore measuring the velocity of the 
host galaxy of one SNe can be as accurate as the veloc- 
ity obtained from 25 galaxies using the D — a relation. 
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FIG. 2: Contribution to the ISW effect power spectra from 
different 3D wavenumbers, k, and thin redshift z shells, for 
various multipoles t. As in Fig. 1, top panels display the 
Tv kernel, amplified by a factor of 100, whereas the bottom 
panels show the TS kernel. Kernels are given by £(£ + 1)/2tt 
times the integrand in eq. (1). 

In turn, the redshift for nearby galaxies (SNe host or 
not) is currently measured using narrow lines from the 
host galaxies with an uncertainty aicz) = 30km. s -1 that 
is dominated by statistical errors. As such for a Hubble 
constant of h = 0.72, errors in redshift measurements will 
be subdominant as compared to errors in Hod for z >0.02 
(0.004) if we use SNe (galaxy) based distance estimators, 
respectively. Since galaxy peculiar radial velocities typ- 
ically have a rms of 300 km.s -1 we reach a signal-to- 
noise ratio S/N ~ 1 per source at z ~ 0.02(0.004). If 
we were to volume-average the distance measurements 
coming from SNe observations over samples at fixed z, 
the systematic error limit would be <j(m) — 0.02 [27]. 
Assuming that we had enough SNe to reach this limit 
within boxes of depth a(cz) — 30km. s" 1 , then one could 
get a 1% distance measurement and the S/N per source 
would reach 1 at z ^ 0.1. 

On the other hand, the kSZ determined peculiar ve- 
locities of galaxy clusters do not rely on a distance 
determination, and associated z-independent errors of 
~ 100 km s -1 may be achievable [28, 29]. However the 
signal is limited to z > 0.1 in order for the primordial 
CMB not to be a dominant source of confusion. Its power 
as a tracer of the large-scale gravitational potential has 
already been studied in [24] and as a DE probe in con- 
junction with other density tracers in [30, 31]. 

The nature of the signal we are investigating demands 
wide but rather shallow surveys in comparison to other 
probes. To illustrate the detectability of this new corre- 
lation, we will consider two ideal full-sky surveys, with a 
distribution of sources with median redshift z m = 0.15 or 
z. m = 0.3, and study their performance in ISW detection 



FIG. 3: Contribution of multipoles to the signal-to-noise ratio 
S/N for Tv (solid) and TS (dashed) for full-sky surveys with 
different depths, for our baseline ACDM model. Total (S/N) 2 
is given by the area under the curves. Velocities give a S/N in 
cross-correlation with the temperature that is systematically 
larger than that for densities, TS, the gain factor depending on 
survey depth and cosmology. Note that S/N scales as yj f s k y , 
where f sky is the fraction of sky covered by the survey. 

and DE constraints as compared to the T5 correlation. 
In particular, we propose to combine full-sky peculiar ve- 
locity surveys from SNe at z < 0.1 and, complementarily, 
kSZ cluster surveys for sources at z > 0.1. We define the 
cross-correlation signal-to- noise ratio S /N as 

(r<TX\2 

{S/N)* = £(2£ + l)f sky {cJx y + ^ TcfX (5) 

where X = v or S. Fig. 3 shows the contribution of differ- 
ent multipoles to the S/N depending on survey depth, for 
a ACDM model with Qde = 0.75 that will be our base- 
line. The area under the curves give the total (S/N) 2 . We 
find that for low z, Tv peaks at similar (albeit slightly 
larger) £'s than TS, i.e, I ~ 10 — 20. For our baseline cos- 
mology one gets, for the z m = 0.3 survey, (S/N)t\, = 5, 
that is 25% larger than the corresponding significance 
for the ISW detection from the T5 correlation. The shal- 
lower survey, z m — 0.15, leads to more moderate sig- 
nificance, (S/N)t„ = 3.2. However, this is about factor 
of 2 larger than that for T5. This comes from the fact 
that the lowest redshift sources give the dominant con- 
tribution to the overall S/N. In other words, as shown by 
Fig. 2, progressively deeper velocity surveys are less op- 
timal for ISW measurements since the Tv signal quickly 
drops with redshift unlike TS. In general, the relative 
gain in S/N when using velocities as compared to densi- 
ties can be understood by computing the ratio of S/N for 
Tv with respect to TS for a given multipole £, that scales 

as (Cj v /Cf 5 )Jcf s /C^ v . Although at low multipoles 
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FIG. 4: Constraints on the DE parameters from the cross- 
correlation of the CMB with shallow full-sky surveys. Con- 
tours show 68% (inner line) and 95% (outer line) confidence 
levels for 1 DOF. Solid (dashed) lines show parameter errors 
from the Tv (T6) correlation, whereas thick (thin) lines cor- 
respond to constraints with (without) PLANCK priors. Left 
(right) panel shows the case for surveys with median red- 
shift z m = 0.15 {z m — 0.3) and top (bottom) panels display 
two different choices of the DE equation of state: w = — 1 
(ACDM) and -2/3 (quintessence). 



the signal is about two orders of magnitude smaller for 
velocities, (Cj v /Cf s ) ~ 0.01 (see Fig.l), the correspond- 
ing ratio of noise terms arising from the auto-correlations 
Cf S l c T k, 1q4 > what makes the relative S/N exceed 
unity. We note that this Tv significance gain over TS 
increases for more strongly DE dominated cosmologies. 
In particular, for a shallow survey (z m = 0.15), velocities 
can measure the cosmic low multipoles with a significance 
> 2 times larger than dark-matter density tracers, such 
as galaxies, in cross-correlation with the CMB. 



COSMOLOGICAL IMPACT AND CONCLUSIONS 

In the context of flat CDM models, an increased signif- 
icance in the ISW measurement from the Tv correlation 
directly translates into tighter DE constraints. Fig. 4 
displays the expected parameter constraints that can 
be achieved for two shallow full-sky surveys. For com- 
parison purposes, we show constraints from Tv (solid 
lines) and the usual TS (dashed) correlations. In gen- 
eral, the Tv estimator yields better DE constraints than 
the TS irrespective of the baseline cosmology used or sur- 
vey depth, provided z m < 0.3, beyond which the usual 
TS cross-correlation outperforms Tv. A shallow velocity 
survey with z m — 0.15 can deliver a (symmetrized) er- 



ror <t(Qde) = 0.065 for our baseline ACDM cosmology, 
that is almost a factor of 2 better accuracy than what 
densities can deliver (see e.g, lower left panel in Fig 4). 
These constraints are comparable to what is expected 
for a(VLr>E) from PLANCK[35] (see comparison between 
thin and thick solid lines in the left panels), although such 
an outstanding CMB dataset will largely help breaking 
the degeneracy with w. For the deeper survey, z rn = 0.3, 
velocities yield a factor of 2 better measurement of the 
DE density (j(Q,de) = 0.03 than the shallower survey 
(solid contours in right panels as compared to same lines 
in left panels), but this is comparable to what a density 
tracer would deliver (dashed contours in corresponding 
panels). Velocities can also help breaking the degener- 
acy with w, but, as expected from a low redshift survey, 
to a much lower extent than PLANCK. Note however 
how the DE errors depend on the baseline cosmology 
used: Tv measurements can provide significanlty sharper 
constraints on quintessence models than on ACDM and 
yield better measurements of £Ide than PLANCK. This 
is shown in the lower panels of Fig 4 where errors from 
the Tv correlation alone (thin solid lines) do not improve 
when adding PLANCK priors (thick solid lines). 

In practice several factors could degrade these opti- 
mal expectations. In cross-correlations S/N cx ag and 
thus current estimates of as — 0.8 would imply a 10% 
weaker detection than we estimated from our baseline 
model. Lensing could introduce noise in Tv as we go 
deeper in redshift although its large scale contribution 
is small. More importantly, we neglected so far the fact 
that we deal with a discrete sampling of the underly- 
ing velocity field. Note that for velocities no shot noise 
affects their measurement, unlike the case for densities. 
However, if we consider the case for a z m = 0.3 survey, 
90% of sources lie below z = 0.8. Since there are around 
10 5 clusters with masses greater than 10 14 h _1 M Q in 
this volume, a full-sky cluster kSZ survey up to z = 0.8 
would provide a good sampling of the field for modes 
k < 0.1 Mpc -1 which is sufficient for the CMB- velocity 
signal according to Fig. 2. 

In conclusion, the use of wide and shallow velocity sur- 
veys can substantially increase our ability to extract the 
ISW signal from a cross-correlation with the CMB and 
thus help other cosmological probes to distinguish grav- 
ity from alternatives. Peculiar velocities are highly cor- 
related with other large-scale structure tracers such as 
the galaxy density distribution but probe larger physi- 
cal scales. In addition, the newly proposed correlation 
is a priori less affected by systematics that plague the 
temperature-galaxy correlation such as galaxy bias and 
shot noise. In particular, combining shallow velocity sur- 
veys with complementary data from planned deep sur- 
veys (such as LSST, DES, Pan-STARRS, WFMOS) can 
fully exploit the ability of tomographic analyses of future 
surveys to reveal the nature of dark-energy and grav- 
ity on the largest physical scales [32, 33]. Furthermore, 
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the temperature-velocity correlation is a more power- 
ful probe of the CMB lowest multipoles than the stan- 
dard correlation with dark-matter density tracers, such 
as galaxies, and hence it could shed new light on inflation 
and the physics of the primordial universe. 

We are grateful to Niayesh Afshordi, Pierre Astier, En- 
rique Gaztanaga, Joe Hennawi, Wayne Hu, Mike Hud- 
son, Arthur Kosowsky, Robert Lupton, Ramon Miquel, 
Jordi Miralda-Escude, Pat McDonald, Bohdan Paczyn- 
ski, David Schlegel, David Spergel, and Martin White 
for useful comments. PF acknowledges support from the 
Spanish MEC through a Ramon y Cajal fellowship, and 
grants AYA2006-06341, 2005SGR-00728. 



* Electronic address: fosalba@ieec.uab.es 
' Electronic address: olivier@cita.utoronto.ca 
[1] R. K. Sachs and A. M. Wolfe, Astrophys. J. 147, 73 
(1967). 

[2] R. G. Crittenden and N. Turok, Physical Review Letters 

76, 575 (1996). 
[3] R. Bean and O. Dore, Phys. Rev. D 69, 083503 (2004). 
[4] A. Lue et al, Phys. Rev. D 69, 044005 (2004). 
[5] J. Garriga et al, Phys. Rev. D 69, 063511 (2004). 
[6] Y.-S. Song, I. Sawicki, and W. Hu, ArXiv Astrophysics 

e-prints (2006), astro-ph/0606286. 
[7] S. Boughn and R. Crittenden, Nature (London) 427, 45 

(2004). 

[8] M. R. Nolta et al, Astrophys. J. 608, 10 (2004). 
[9] P. Fosalba and E. Gaztanaga, MNRAS 350, L37 (2004). 
[10] P. Fosalba et al, Astrophys. J. Lett. 597, L89 (2003). 
[11] R. Scranton et al, ArXiv Astrophysics e-prints (2003), 

astro-ph/0307335. 
[12] N. Afshordi, Y.-S. Loh, and M. A. Strauss, Phys. Rev. D 

69, 083524 (2004). 
[13] A. Cabre et al, ArXiv Astrophysics e-prints (2006), 
astro-ph/0603690. 



[14] N. Afshordi, Phys. Rev. D 70, 083536 (2004). 

[15] A. Cabre, P. Fosalba, E. Gaztanaga, and M. Manera, 

ArXiv Astrophysics e-prints (2007), astro-ph/0701393. 
[16] M. A. Strauss and J. A. Willick, Phys. Rep. 261, 271 

(1995). 

[17] M. J. Hudson et al, in ASP Con]. Ser. 201: Cosmic 

Flows Workshop, edited by S. Courteau and J. Willick 

(2000), pp. 159-+. 
[18] D. Sarkar et al, ArXiv Astrophysics e-prints (2006), 

astro-ph/0607426. 
[19] A. G. Riess et al, Astrophys. J. Lett. 488, L1+ (1997). 
[20] P. Astier et al, A&A 447, 31 (2006). 
[21] R. A. Sunyaev and I. B. Zeldovich, MNRAS 190, 413 

(1980). 

[22] T. Haugboelle et al, ArXiv Astrophysics e-prints (2006), 

astro-ph/0612137. 
[23] B. A. Benson et al, Astrophys. J. 592, 674 (2003). 
[24] O. Dore et al, Astrophys. J. Lett. 585, L81 (2003). 
[25] E. T. Vishniac, Astrophys. J. 322, 597 (1987). 
[26] N. Kaiser, Astrophys. J. 282, 374 (1984). 
[27] A. G. Kim et al, MNRAS 347, 909 (2004). 
[28] D. Nagai, A. V. Kravtsov, and A. Kosowsky, Astrophys. 

J. 587, 524 (2003). 
[29] G. P. Holder, Astrophys. J. 602, 18 (2004). 
[30] O. Dore, J. F. Hennawi, and D. N. Spergel, Astrophys. 

J. 606, 46 (2004). 
[31] S. DeDeo et al, ArXiv Astrophysics e-prints (2005), 

astro-ph/0511060. 
[32] W. Hu, Phys. Rev. D 65, 023003 (2002), astro- 

ph/0108090. 

[33] M. Kesden, M. Kamionkowski, and A. Cooray, Physical 
Review Letters 91, 221302 (2003). 

[34] we ignored here the velocity of the observer since it is 
well measured through the CMB dipole 

[35] we use Fisher matrix priors from PLANCK temperature 
and polarization data, as given by [32], i.e, we assume 
o{Vt DE ) = 0.1, a(w) = 0.32 for ACDM with {l DE = 0.75, 
and o{Q,de) = 0.2, a(w) = 0.5 for a quintessence model 
with w — —2/3 and Q.de = 0.75. 



